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ABSTRACT: The characteristics of whole waveforms of 
magnetizing inrush current and fault differential current were 
analyzed. By virtue of the similarity difference between the two 
currents and combining the superiorities of the Hausdorff 
distance algorithm, a novel criterion of the transformer 
differential protection based on the Hausdorff distance 
algorithm was proposed. The Hausdorff distance between the 
normalized differential current sequence and standard 
sinusoidal wave sequence was applied to distinguish between 
inrush current(including symmetrical inrush current) and fault 
differential «current(including fault differential current 
superimposed with inrush current). The effectiveness, rapidity 
and anti-interference performance of the proposed criterion 


were validated with simulation tests. 
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摘要 : 对 变压器 励磁 涌流 和 故障 电流 波形 整体 形态 的 特征 
进行 分 析 ， 利 用 两 者 之 间 相 似 度 差异 ， 结 合 Hausdor 任 距离 
算法 在 图 形 相似 度 方面 的 优势 ， 构 造 基于 Hausdor 任 距离 算 
法 的 变压器 差 动 保 护 新 判 据 ， 该 判 据 采 用 标 么 化 后 的 差 流 
序列 和 标准 正弦 波 序列 之 间 的 Hausdorff 距离 值 进行 判别 ， 

以 区 分 变压器 的 励磁 涌流 (包括 对 称 性 涌流 ) 和 故障 差 流 
(包括 故障 电流 蕉 加 励磁 涌流 ) 。 通 过 仿真 试验 验证 所 提 新 
判 据 的 有 效 性 、 快 速 性 和 抗 干扰 性 。 


关键 词 : 变压器 ， 差 动 保护 ， 对 称 性 涌流 ， 图 形 相似 度 
Hausdor 任 距离 算法 
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1 引言 

随 着 特 高 压 输电 技术 的 发 展 和 应 用 ， 更 高 电 
压 等 级 输电 工程 的 接 入 使 得 变压器 传输 容量 巨大 
例如 ，2014 年 1 月 27 日 投 运 的 哈密 南 -郑州 
+800kV 特 高 压 直 流 输电 工程 ， 其 额定 输送 功率 达 
8000MW， 大 型 变压器 保护 的 误 动 或 拒 动 都 可 能 带 
来 无 法 承受 的 后 果 ， 因 此 变压器 所 配置 保护 的 性 
能 至 关 重 要 。 差 动 保护 由 于 原理 简单 、 动 作 灵敏 可 
靠 ， 被 广泛 应 用 于 作为 变压器 的 主 保护 。 但 是 ， 由 
于 变压器 铁 芯 存在 饱和 现象 ， 变 压 器 经 历 空 载 合 
闻 和 外 部 故障 切除 时 会 产生 励磁 涌流 ， 这 是 影响 
变压器 差 动 保护 动作 可 靠 性 的 最 主要 的 问题 。 长 期 
以 来 ， 研 究 人 员 一 直 致 力 于 研究 更 优化 的 励磁 涌 
流 识别 方法 ， 并 应 用 于 实际 保护 系统 中 。 
目前 ， 励 磁 涌流 的 识别 方法 主要 分 为 两 大 类 
一 类 是 以 励磁 涌流 的 产生 机 理 为 出 发 点 ， 基 于 理 
论 分 析 提 出 的 识别 方案 ， 如 磁 通 特性 法 "33、 等 效 瞬 
时 电感 /动态 阻抗 法 B4 等 原理 ， 来 提升 差 动 保护 对 
励磁 涌流 的 识别 能 力 ， 但 此 类 方法 必须 用 到 电压 
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测量 量 ， 增 加 了 电压 互感 器 〈PT) 的 投入 成 本 ， 
保护 配置 复杂 化 的 同时 也 带 来 了 PT 断 线 引起 的 保 
护 不 正确 动作 的 风险 ， 实 用 性 受到 限制 。 
男 一 类 则 是 单纯 以 励磁 涌流 的 波形 特征 为 出 
发 点 ， 通 过 波形 的 对 比分 析 来 实现 涌流 的 识别 ， 
例如 基于 二 次 谐 波 特性 已、 波形 对 称 性 四 和 间断 角 特 
性 中 等 判 据 ， 在 现场 得 到 广泛 应 用 。 但 是 由 于 现代 


We 


变压器 铁心 材料 的 改进 、 系 统 中 并 联 电容 器 的 应 用 ， 


以 及 交 直 流 混 联 耦 合 使 得 变压器 工作 的 电磁 暂 态 
环境 更 为 复杂 ， 可 能 使 涌流 和 故障 电流 的 波形 及 
谐 波 特征 差异 弱化 ， 因 涌流 对 称 性 好 而 导致 二 次 
谐 波 制 动 判 据 失 效 差 动 保护 发 生 误 动 [8]， 甚 至 出 
现 现场 二 次 谐 波 制 动 原理 和 波形 对 称 性 制 动 原理 
的 两 套 保护 同时 误 动 的 案例 中。 

值得 注意 的 是 ， 从 波形 特征 来 讲 ， 涌 流 的 产 
生 受 变压器 铁 芯 饱和 的 影响 ， 波 形 与 主要 呈正 弦 


征 点 ， 如 Hausdor 任 距离 算法 上 19。 
对 于 两 组 点 集 4={a.….,ay}，B={b1...,by}， 
Hausdorff 距离 定义 为 : 
H(AB)=max(h(AB)NBA) Ul) 


un 


其 中 : 
h(4,B)= max minla.-b, O) 
h(B,A)= max min Db -a (3) 


bEB aEA 


| 是 点 集 4 和 点 集 B 间 的 距离 范 数 。 式 (2) 中 首 


先 对 点 集 4 中 每 个 点 (如 aj) 与 点 集 B 中 所 有 点 
进行 距离 比较 ， 找到 对 于 a 最 近 的 点 b: 


a —b; fla, b | (1<k<9 [1)，j|a 一 即 为 对 
应 a; 点 的 最 小 距离 ，h(4,B) 即 为 所 有 和 集合 4 中 点 的 
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特性 的 故障 电流 存在 差异 ， 因 此 可 利用 这 样 的 差 
异 来 构造 波形 相似 度 的 判 据 ， 此 类 判 据 简单 易 实 
现 ， 例 如 ， 文 献 [10-14] 采 用 不 同 算法 来 对 差 流 波 
形 正 弦 度 或 波形 相关 系数 进行 计算 ,设置 门槛 ， 
以 有 效 区 分 涌流 和 故障 电流 。 但 上 述 算法 均 需 首先 
对 差 流 直流 分 量 进 行 滤波 处 理 ， 其 次 ， 在 执行 计 
算 时 ， 数 据 窗 需要 至 少 半 个 甚至 一 个 周波 才能 做 
出 判断 ; 并 且 ， 在 差 流 采样 序列 中 有 幅 值 较 大 的 
畸变 点 时 ， 计 算 结 果 将 受到 很 大 影响 ， 算 法 抗 干 
扰 能 力 较 弱 。 因 此 ， 我 们 考虑 进 
判别 快速 并 具有 较 强 抗 干扰 性 的 波形 相似 度 识 别 
新 判 据 ， 以 提升 变压器 差 动 保护 的 动作 性 能 。 

本 文 在 对 变压器 励磁 涌流 〈 包 括 对 称 性 涌 
流 ) 与 故障 电流 〈 包 括 故障 电流 县 加 典型 涌流 ) 
波形 差异 分 析 的 基础 上 上， 结合 Hausdorff 距离 算法 
在 时 窗 选 择 、 频 率 适 应 性 及 抗 数 据 丢 失 能 力 等 方 盏 
的 优势 ， 构 造 基于 Hausdorff 距离 算法 的 变压器 差 
动 保护 新 判 据 ， 对 变压器 经 历 空 载 合 闸 、 带 故障 合 
闻 、 区 内 故障 和 涌流 后 又 发 生 故 障 等 各 种 情况 时 ， 
新 判 据 的 有 效 性 进行 了 仿真 试验 验证 ， 结 果 表 明 
该 判 据 原理 简单 ， 易 于 实现 ， 能 对 变压器 励磁 涌 
流 和 故障 差 流 进行 正确 识别 ， 在 速 动 性 和 抗 干 扰 
性 方面 优势 明显 。 


2 Hausdorff 距离 算法 的 基本 原理 


图 形 相似 度 算法 是 一 种 比较 两 幅 图 像 的 整体 
特性 而 不 局 限于 某 些 小 范围 特征 影响 的 算法 。 图 形 


步 研 究 算 法 简单 、 


最 小 距离 的 最 大 值 ， 称 为 从 点 集 4 到 点 集 B 的 单 
向 Hausdorff 距 离 。 式 (1) 的 Hausdor 企 距离 为 单 向 距 
离 h(4,B) 和 有 (B,4) 的 较 大 值 。 

Hausdorff 距离 是 描述 两 组 点 集 之 间 相 似 程度 
的 一 种 量度 ， 它 度量 了 两 组 点 集 间 的 最 大 不 匹配 
程度 ， 亦 即 反 映 目标 图 形 与 模板 图 形 边缘 特征 点 
之 间 的 匹配 程度 。 在 生物 医学 领域 ， 己 有 研究 采用 
Hausdorff 距离 算法 来 实现 生物 医学 信号 的 相似 性 
判断 ， 例 如 利用 Hausdorff 距离 算法 来 实现 对 随时 
间 变 化 的 心 电 波形 与 正常 心 电 波形 相似 度 的 判别 


[17,18] 


一 让 


而 在 继 电 保 护 系 统 中 ， 传 感 器 所 采集 的 电流 
数据 以 二 维 点 集 的 形式 存在 ， 可 看 作 是 一 个 以 时 
间 为 横 坐 标 ， 幅 值 大 小 为 纵 坐 标的 离散 时 间 序 列 ， 
每 一 个 电流 数据 点 都 相当 于 图 形 的 某 个 特征 点 ， 
这 与 前 面 所 提 心 电波 形 类 似 ， 因 此 也 可 以 直接 用 
于 进行 Hausdorff 距离 计算 。 

3 Hausdorff 距离 算法 的 特点 及 对 差 动 保护 需 
求 的 适应 性 

相 比 于 传统 差 动 保护 算法 中 应 用 的 逐 点 比较 

差 值 寻 获 的 算法 ，Hausdorff 距离 算法 具有 以 下 特 
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第 一 ，Hausdorff 距离 算法 可 以 根据 需求 自由 
选择 数据 窗 长 度 。 传 统 的 保护 算法 一 般 采 用 健 里 叶 
算法 ， 常 用 的 有 全 周 健 里 叶 算 法 以 及 半 周 傅 里 叶 
算法 。 但 是 这 两 种 算法 都 需要 保证 时 间 窗 的 长 度 是 


识别 算法 一 般 分 为 两 大 类 ， 一 类 基于 图 形 灰 度 信 
息 ， 如 PCA 主 元 分 析 算 法 ;， 另 一 类 基于 图 形 的 特 


电流 半 周 期 的 整数 倍 。 在 继 电 保护 工程 应 用 实际 中 ， 
这 种 时 间 窗 特性 导致 处 理 故 障 信号 时 至 少 有 10ms 


的 延迟 ， 并 且 不 能 根据 需求 灵活 设置 时 间 窗 长 度 。 
而 Hausdorff 距离 算法 不 涉及 信号 从 时 域 到 频 域 的 
投射 ， 因 此 时 间 窗 的 设置 可 以 更 灵活 ， 考 虑 到 继 
电 保护 的 速 动 性 需求 ， 在 下 面 的 算法 设计 时 ， 拟 
采用 1/4 周波 的 Hausdorff 距离 计算 时 窗 ， 其 窗 长 
仅 为 DFT 的 1/4 到 1/2， 进 一 步 体现 其 对 继 电 保护 
速 动 性 需求 的 适应 性 ; 

第 二 ，Hausdorff 距 离 算 法 受 采 样 频率 的 影响 
较 小 。 传 统 傅 里 时 算法 对 保护 装置 采样 频率 的 要 求 
较 高 ， 如 果 保 护 装 置 采样 率 过 低 ， 会 导致 信号 在 
频 域 的 投射 不 准确 。 而 Hausdorff 距离 算法 多 考虑 
的 是 整体 特征 的 一 致 性 ， 而 对 采样 点 的 时 域 对 齐 
及 采样 率 的 统一 性 并 没有 严格 要 求 ， 例 如 1.2kHz 
的 采样 率 相 比 于 4kHz 的 采样 率 ， 只 是 波形 特征 值 


因为 比较 的 是 差 流 波形 形态 的 特征 ， 其 幅 值 
的 影响 在 Hausdor 作 距离 计 算 中 应 当 被 剔除 .剔除 的 
方式 为 : 在 一 定 的 采样 率 下 ， 首 先 利 用 一 个 很 短 
的 时 间 窗 获取 差 流 的 幅 值 信 息 ， 之 后 以 该 幅 值 信 
息 作为 基准 ， 对 差 流 进行 比例 压缩 和 纵向 平移 ， 

得 到 标 和 后 的 差 流 序列 。 标 么 后 的 差 流 序列 波形 幅 
值 变化 范围 设 定 为 [-1,1]。 将 标 乏 后 的 差 流 序列 作为 
Hausdorff 距离 算法 目标 图 形 的 边缘 特征 点 ， 将 相 
同 采样 频率 的 幅 值 为 1 的 标准 正弦 波 序列 作为 
Hausdorff 距离 算法 模板 图 形 的 边缘 特征 点 ， 计 算 
两 者 间 的 Hausdor 人 ff 距离。 波形 经 过 标 乏 化 处 理 后 ， 

该 Hausdorff 距离 计算 值 必然 落 在 [0,1] 之 间 ， 数 值 
越 小 ， 代 表 差 流 序列 的 波形 越 接近 正弦 波 ， 反 之 ， 
数值 越 大 ， 表 示 差 流 序列 与 正弦 波 的 相似 度 越 低 。 


的 选取 比较 稀疏 ， 不 会 造成 特征 值 的 错误 ， 对 于 
波形 整体 特征 的 判断 没有 影响 ; 

第 三 ，Hausdorff 距离 算法 具有 抗 数 据 丢 失 能 
力 。 个 别 数据 点 的 丢失 往往 会 造成 采样 值 差 动 算法 
失效 。 但 对 于 Hausdorff 距离 算 法 而 言 ， 并 不 影响 
其 对 图 形 整 体 特征 的 判断 ， 利 用 该 特性 ， 在 进行 
距离 计算 前 ， 可 根据 实际 情况 设计 舍弃 若干 极 值 
点 ， 并 不 影响 整体 距离 值 计 算 结 果 ， 对 于 差 流 序 
列 中 的 随机 噪声 具有 很 强 的 抗 干扰 性 。 


4 基于 Hausdorff 距 离 算法 的 差 流 波形 相似 性 
判断 


里 论 上 ， 对 于 理想 故障 电流 而 言 ， 不 考虑 非 
周期 分 量 及 其 幅 值 变化 ， 其 波形 基本 呈现 正 弱 波 
特征 ， 而 对 于 涌流 来 讲 《〈 单 向 或 对 称 性 ) ， 因 其 
产生 受 变 压 器 铁 芯 饱和 的 影响 ， 幅 值 的 上 升 存在 
一 个 加 速 ， 所 以 呈现 尖 波 的 形态 ， 与 正弦 波 存在 
很 大 差异 。 因 此 我 们 可 以 将 正弦 波 作 为 基准 波形 ， 
将 采样 得 到 的 差 流 波形 与 基准 波形 进行 相似 度 的 
判断 ， 帮 接近 基准 正弦 波 ， 则 判别 为 故障 差 流 ， 
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若 偏离 基准 正弦 波 超过 一 定 程度 ， 则 认为 是 涌流 ， 


以 此 来 决定 闭锁 还 是 开放 差 动 保护 。 

基于 上 述 讨 论 ， 可 以 设计 这 样 的 比较 过 程 : 
将 目标 波形 ( 差 流 ) 与 一 个 按照 某 种 方式 构造 出 
来 的 标准 正弦 波 基准 波形 (模板 〉 进行 相 似 度 或 
匹配 度 的 评估 ， 来 决定 其 所 对 应 的 扰动 到 底 是 内 


里 论 上 ， 阁 是 内 部 故障 引起 的 差 流 ， 标 乏 后 的 故 
障 差 流 序列 与 幅 值 为 1 的 标准 正弦 波 序列 的 
Hausdorff 距离 应 接近 0; 而 对 于 涌流 引起 的 差 流 
情况 ， 两 者 Hausdorff 距 离 应 较 大 ,我 们 首先 考察 几 
组 理想 故障 差 流 和 涌流 与 标准 正弦 波 序列 的 
Hausdorff 距离 计算 值 ， 再 根据 具体 情况 制定 合理 
的 整定 原则 ， 来 对 故障 和 涌流 的 情况 进行 有 效 区 
pa 
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1 标 么 化 差 流 与 标准 正弦 波 相 似 性 比较 
Fig.1 Comparison of normalized differential currents with 
standard sinusoidal waveform 
图 1(a)~(d) 分 别 给 出 了 标 么 化 处 理 后 的 变压器 
内 部 故障 差 流 、 空 载 合 闸 单 向 励磁 涌流 、 对 称 性 涌 
流 、 带 故障 合 曾 差 流 故障 电流 县 加 涌流 ) 与 标准 


部 故障 还 是 涌流 。 根据 前 文 介 绍 ，Hausdor 任 距离 算 
法 是 反映 目标 图 形 与 模板 图 形 边缘 特征 点 之 间 匹 
配 程度 的 算法 ， 因 此 ， 应 用 其 进行 差 流 波 形 与 模 
板 间 的 相似 性 判断 。 下 面 对 基 于 Hausdorff 距离 算 
法 的 波形 相似 性 判断 逻辑 进行 设计 。 


正弦 波 的 相似 性 比较 (图 中 虚线 为 标准 正弦 波 ， 
实 线 为 差 流 波形 ) ， 差 流 序列 采样 频率 为 4kHz， 
即 每 周波 采样 80 个 点 。 对 上 述 四 组 波形 进行 
Hausdorff 距离 计算 ，Hausdorff 距离 算法 时 间 窗 取 
为 /4 个 周波 ， 距 离 计 算 值 用 刀 表 示 。 表 1 分 别 列 


举 了 一 
个 周波 结束 时 计算 


8 的 均值 。 


需要 指出 的 是 ， 以 1/4 周波 滑动 时 
1(b) 的 空 龙 
流 ， 在 第 3/4 周波 和 第 4/4 周波 期 间 ， 
无 法 形成 标准 正 
这 种 情况 ， 定 义 五 值 为 1。 


电流 采样 数据 时 ， 对 于 对 应 于 图 
向 励磁 涌 
为 时 间 窗 内 的 数据 均 接近 于 零 ， 
弦 模 板 的 幅 值 ， 对 于 


> 


个 完整 周波 内 ， 第 1/4、 第 1/2、 第 3/4 和 第 1 


间 窗 去 截取 
合 闸 单 


表 1 标 么 化 差 流 与 标准 正弦 波 Hausdorff 距 离 计算 值 


Tab.1 Hausdorff distance values between the normalized 


differential current with 人 sinusoidal waveform 


差 流 类 型 
EE 14 | 第 12 | 第 3/4 | 第 1 
周波 | 周波 | 周波 | 周波 


内 部 故障 单纯 故障 差 流 (a) 
磁 涌 流 (b) 


0.482 0.467 


0.097 0.041 | 0.082 | 0.068 
0.582 0.563 | 0.580 | 0.544 


0.118 | 0.112 


从 表 1 可 以 看 到 : 


中 对 于 变压器 内 部 故障 ， 标 么 化 差 流 序列 与 
标准 正弦 波 相似 性 非常 好 ，4 个 1/4 周波 窗 计算 的 


五 值 非常 低 ， 接 近 0; 


包 对 于 单 向 励磁 涌流 和 对 称 性 ; 
较 高 ， 这 表明 两 者 与 


算出 的 瑟 值 均 超过 0.5， 数 值 
标准 正弦 波 的 相似 度 低 。 


故障 电流 与 涌流 的 车 加 ， 如 图 
与 标准 正弦 波 的 相似 性 
在 0.4-0.5 之 间 ， 但 在 涌流 幅 
流 已 开始 呈现 典型 故障 电流 特 4 
1 个 周波 内 ， 差 流 
五 值 在 0.12 以 下 。 
基于 上 述 分 析 ， 
为 判 据 整 定 的 依据 ， 但 
差 流 与 标准 正弦 波 之 间 


将 内 部 故障 时 的 理 
因 理 想 情况 下 内 部 故障 时 ， 
五 值 理论 上 接近 0 


涌 流 序列 ， 讨 : 


@ 对 于 变压器 带 故 障 合 闻 的 情况 ， 其 差 流 为 
1(d) 所 示 ， 差 流 序列 
在 前 1/2 周波 内 较 大 ， 
直 较 小 的 负 半 波 ， 差 
正 ， 即 在 第 
阐 列 与 标准 正弦 波 的 相似 度 增 大 ， 


瓦 值 


3/4- 到 第 


论 五 值 作 


甚至 即 


为 0， 工 程 实际 中 ， 若 以 0 


值 作 为 整定 基准 ， 则 无 


法 有 效 进行 可 靠 系数 的 乘除 和 灵敏 度 校 验 。 考 虑 到 


标 乏 化 处 理 后 ,五 计算 值 必然 落 在 [0,1] 之 间 ， 
作为 判 据 基准 ， 设 定 
.， 即 在 计算 标 终 化 故障 差 


我 们 采用 五 计算 值 的 补 重 
HS(D=1-H(O, =1,2,3.. 


流 序 列 与 幅 值 为 1 的 标准 正弦 波 序列 的 Hausdorff 
距离 值 厂 后， 再 计算 出 对 应 的 值 ， 


因此 


以 HS 值 大 


小 作为 判别 内 部 故障 和 励磁 涌流 的 依据 。 


因此 ， 整 定 原则 如 下 : 
Lu Heory /K, rel 
其 中 ，K 一般 取 1.15~1.3 之 间 。 


=]， 不 妨 取 开 Rer1.3， 则 到 Ss =0.77。 


别 结果 ， 如 表 2 所 示 。 
由 表 2 可 知 ， 对 于 内 部 故障 ，H5 在 整个 时 间 


窗 扫 
流 ，HS 均 小 0.5。 对 于 带 故 障 合 
前 半 个 周波 HS 均 小 于 动作 门槛 ,1 
已 经 提升 到 门槛 之 上 ， 达 到 了 0.88 左右 ， 


于 门槛 值 。 


根据 表 1 的 数据 计算 出 相应 的 到 值 


(5) 
由 于 HS ineory 


和 判 据 关 


的 范围 内 均 大 于 0.9， 而 对 于 两 类 典型 的 涌 
闸 的 案例 ， 虽 然 在 


在 后 


二 下 


王者 局 


半 周 波 
» 
三 | 

NT 


表 2 标 么 化 差 流 与 标准 正弦 波 的 HS 计算 值 和 判别 结 


Tab.2 HS values and discrimination results between the 


normalized differential current with standard sinusoidal 


还 是 对 称 性 ; 
保护 能 可 靠 闭 锁 ; 


内 部 故 
障 差 流 (a) 


空 载 合 闸 单 向 励 


磁 涌 流 (b) 


涌流 (c) 
故障 电 
流 (d) 


差 流 类 型 
障 征 EE 
闸 对 称 性 


流 全 加 涌 


0 


| 5 值 (457077) | (HS;e=0.77) 


第 4 | 第 12 | 第 34 | 第 1 
周波 | 周波 | 周波 | 周波 


0.903 0.951 0.918 | 0.932 


>0.77， 保 护 动 作 
| 


<0.77， 保 护 闭锁 


0.418 0.437 0.420 | 0.456 


<0.77， 保 护 闭锁 


0.518 0.533 0.882 | 0.888 


<0.77， 


保护 闭锁 | >0.77， 


保护 动 
作 


综 上 可 知 ， 对 于 涌流 情况 ， 无 论 


对 于 


要 1/4 个 周波 ， 即 Sms， 


使 保护 正确 动作 ;对 于 
最 迟 3/4 个 周 


本 判 据 ， 
判断 ， 开 放 保护 。 


是 单 向 涌流 


涌流 ，H5S 值 始终 在 0.77 的 门槛 值 以 下 ， 
内 部 单纯 故障 
该 判 据 即 可 做 出 正确 判断 ， 
故障 差 流 竺 


差 流 ， 只 需 


涌流 ， 采 用 


加 


波 ， 即 15ms 也 可 


做 出 


E 确 


5 基于 Hausdorff 距离 算法 的 变压器 差 动 保护 
新 判 据 


别 变 压 器 的 内 部 故障 
基于 


序列 的 整体 特 行 
域 到 频 域 的 投射 ， 


时 间 窗 的 选取 灵活 


根据 前 面 的 分 析 ，Hausdorff 距离 算法 可 以 判 
差 流 和 励磁 涌流 。 
FE 进行 判断 ， 不 涉及 信号 从 时 


该 算法 主要 


闲 


在 


于 差 


动 保 护 时 ， 为 保证 获取 到 周期 性 序列 的 极 值 点 ， 
用 于 序列 波形 标 和 化 处 理 ， 可 将 时 间 窗 设置 为 1/4 
周波 。 从 第 1/4 个 周波 开始 ， 时 间 窗 每 向 后 移动 一 
个 采样 点 ， 都 可 以 更 新 一 个 Hausdorff 距 离 值 ， 实 
际 将 生成 一 个 五 值 序列 ， 进 而 产生 HS 值 序列 ， 用 
于 实时 判断 差 流 变化 暂 态 特征 。 

根据 上 一 小 节 的 分 析 ， 对 于 一 般 性 的 变压器 
内 部 故障 差 流 或 励磁 涌流 〈 包 括 对 称 性 涌流 ) ， 
该 算法 在 1/4 周波 〈$ms) 即 可 做 出 正确 判断 ， 而 
对 于 特殊 场景 ， 如 故障 电流 登 加 清流 的 情况 ， 最 


一 让 


述 3/4 周波 〈1$ms) 便 可 做 出 响应 使 保护 正确 动作 。 


这 体现 了 该 算法 的 快速 性 。 据 此 构造 基于 Hausdorff 
距离 算法 的 变压器 差 动 保护 新 判 据 ， 流 程 如 图 2 所 


a 
N 一 莹 流 幅 值 、、 
是 否 越 限 一 
折 
启动 判 据 


机 


提取 1/4 周 波 数据 窗 


获取 差 流 序列 幅 值 


SS 
一 柱 信 > 门槛 信 
3 
标 么 化 差 流 序列 
eg 励磁 涌流 
计算 标 乏 化 差 流 序列 [| 闭锁 保护 
与 标准 正弦 波 的 豆 值 
v 


计算 Hs 值 


IN 


Ss 
—H5 >HSsef? 一 | 


Y 
YY 
内 部 故障 ， 保 护 动作 


2 基于 Hausdor 企 距离 算法 的 差 动 保护 判 据 流程 图 


Fig.2 Flow chart of the transformer differential protection 


criterion based on the Hausdorff distance algorithm 


需要 指出 的 是 ， 该 判 据 是 基于 差 动 保护 差 流 
值 越 限 启动 后 才 进 行 判别 ， 现 场 差 流 幅 值 的 获 
目前 至 少 半 个 周波 ， 而 这 半 个 周波 的 数据 同时 
以 拿 来 用 作 该 判 据 第 一 个 幅 值 的 获取 ， 在 第 一 
昼 值 获取 后 ， 我 们 就 可 以 往 前 追溯 相应 采样 点 
同步 生成 标准 正弦 波 序 列 ， 来 对 标 乏 化 差 流 波 
形 进 行 对 比 ， 随 着 时 间 窗 推移 不 停 的 有 采样 点 退 
出 时 间 窗 而 新 的 采样 点 进入 时 间 窗 ， 会 实时 地 将 
新 采样 点 纳入 幅 值 计算 并 形成 新 的 基准 进行 标 乏 
化 处 理 ， 理 论 上 ， 该 判 据 最 快 的 动作 时 间 能 达到 
启动 判 据 同 样 的 动作 速度 ， 并 且 即 使 不 能 立刻 判 


天 祝 到 淄 剖 


断 ， 采 用 1/4 周波 时 间 窗 往 后 推移 ， 加 之 前 面 至 少 
半 个 周波 差 流 幅 值 越 限 启动 判 据 已 存在 的 采样 点 ， 
也 是 够 做 出 判别 ， 一 般 情况 下 ， 在 判 据 启 动 后 5ms 
能 做 出 判断 。 

在 计算 量 方面 ， 以 4kHz 采样 率 为 例 ，1/4 周 
波 窗 长 含 20 个 采样 点 ， 完 成 一 次 Hausdor 人 f 距 离 计 
算 总 共 执 行 约 2000 次 加 减 运算 ， 这 对 于 目前 功能 
强大 的 微机 来 说 ， 所 耗 时 间 微 乎 其 微 ， 不 会 对 差 
动 保护 的 判别 造成 附加 的 延 时 。 

综 上 ， 该 判 据 具 有 很 高 的 动作 可 靠 性 ， 因 此 ， 
在 用 于 大 型 变压器 差 动 保护 时 可 采用 分 相 制 动 的 
方式 ， 与 传统 二 次 谐 波 制 动 判 据 一 起 ， 构 成 原理 
不 同 的 两 套 主 保护 ， 以 满足 保护 双重 化 配置 的 要 
求 。 


6 判 据 仿真 试验 验证 


本 小 节 对 上 述 基 于 Hausdorff 距离 算法 的 变 压 
器 差 动 保护 新 判 据 进行 仿真 试验 验证 。 采用 文献 [8] 
所 述 +800kV 特 高 压 直 流 输电 工程 仿真 模型 ， 以 整 
流 侧 换 流 站 极 I 高 端 一 组 变压器 为 研究 对 象 ， 将 变 
压 器 一 、 二 次 侧 电流 全 部 折算 到 二 次 侧 形 成 各 种 扰 
动 场景 下 该 组 大 差 保 护 差 流 ， 仿 真 时 长 为 1s， 差 
流 和 标准 正弦 波 序列 采样 频率 为 4kHz， 即 每 周波 
采样 80 个 点 ; 差 流 幅 值 越 限 门槛 采用 常规 的 
0.25p.u.，Hausdorff 距离 算法 时 间 窗 取 为 1/4 周波 ， 
动作 门槛 HS,=0.77。 

下 文 给 出 几 例 典型 扰动 的 算 例 仿真 结果 ， 包 
括 变压器 经 历 内 部 故障 、 空 载 合 闸 单 向 励磁 涌流 、 
空 载 合 闸 对 称 性 涌流 、 带 故障 合 闸 、 正 常 空 载 合 闸 
后 发 生 内 部 故障 以 及 变压器 区 外 故障 转 区 内 故障 
的 情况 ， 并 对 算法 抗 干扰 性 进行 了 算 例 验证 。 在 进 
行 Hausdor 任 距离 判 据 HS 值 计 算 时 ， 涌 流 中 对 应 
于 间断 角 部 分 的 点 ， 根 据 上 述 的 处 理 原则 ， 直 接 
赋 五 值 为 1， 对 应 的 到 值 为 0。 
算 例 1: 变压器 正常 内 部 故障 场景 ， 
时 Y/Y 变压器 一 次 侧 出 口 三 相 接 地 故障 。 
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(0) 判 据 到 S 值 计算 序列 
3 差 流 波形 和 判 据 HS 值 计 算 结 果 〈 算 例 1) 


Fig.3 Differential current and value of HS for the proposed 


criterion(Scenario 1) 


如 图 3 所 示 ， 图 (a) 为 差 流 序列 ij， 图 (b) 给 出 
了 仿真 时 间 0.3~0.6s 之 间 标 夭 化 处 型 
zoom 《〈 实 线 表示 ) 和 标准 正弦 波 序列 i 点 状 线 表 
示 ) 的 比较 ; 图 (c) 则 是 根据 图 (b) 所 示 两 个 序列 采 
用 Hausdorff 距离 算法 判 据 后 计算 出 的 到 $ 值 序列 。 
， 在 t=0.41s， 即 故障 发 生 
约 1/4 个 周波 ， HS 计算 值 便 超 过 动作 门槛 
0.77， 保 护 迅 速 正确 动作 。 
算 例 2: 该 算 例 为 单 向 励磁 涌流 场景 ， {=0.4s 
合 闸 空 载 ， 即 A 相合 闸 初 相 角 为 -120"， 合 闸 前 
Y/A 和 Y/Y 变 各 相 初 始 剩 磁 均 为 0。 
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4 差 流 波形 和 判 据 HS 值 计 算 结 果 〈 算 例 2) 


criterion(Scenario 2) 


Fig.4 Differential current and value of HS for the proposed 


算 例 3: 该 算 例 为 对 称 性 涌流 场景 ，t=0.4113s 
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(d) 差 流 二 次 谐 波 百分比 分 析 
5 差 流 波形 和 判 据 HS 值 计算 结果 〈 算 例 3) 


= 


criterion(Scenario 3) 


即 A 相合 闻 初 相 角 为 -120*， 合 闻 前 
Y/A 和 Y/Y 变 A 相 剩 磁 分 别 为 0.85p.u. 和 -0.85p.u.， 
其 他 相 剩 磁 均 为 0。 


Fig.S Differential current and value of HS for the proposed 


算 例 2 和 算 例 3 的 仿真 结果 分 别 如 图 4 和 图 5 
所 示 ， 图 中 各 波形 和 变量 含义 与 图 3 相同 ， 不 效 述 。 
可 以 看 到 ， 无 论 是 算 例 2 的 单 向 励磁 涌流 ， 还 是 算 
例 3 的 对 称 性 涌流 ， 标 乏 化 后 差 流 与 标准 正弦 波 的 
Hausdorff 距离 判 据 HS 值 均 稳定 地 低 于 0.77 的 动 
作 门 槛 值 ， 因 此 判别 为 涌流 情况 ， 能 够 可 靠 闭锁 


EF 


保护 。 值 得 注意 的 是 ， 对 算 例 3 中 差 流 二 次 谐 波 
基 波 百分比 进行 分 析 ， 如 图 5(d) 所 示 ， 可 以 看 到 对 
涌流 的 二 次 谐 波 百 分 比 低 于 通常 设 定 的 15% 
出 动 门槛 ， 关 采用 二 次 谐 波 制 动 判 据 ， 差 动 保 护 
将 误 动 ， 而 采用 本 文 所 提 判 据 则 能 有 效 防止 此 类 


口 


误 动 的 发 生 。 

算 例 4: 该 算 例 为 故障 县 加 涌流 场景 ， 换 流 变 
在 t=0.4s 空投 ， 其 中 Y/Y 换 流 变 带 A、C 两 相 高 阻 
接地 故障 ， 两 台 换 流 变 各 相 初 始 剩 磁 均 为 0， 故障 
接地 电阻 为 70Q。 

4 


00 0 03 04 


03 06 07 08 05 


pe 


0 0 02 03 04 03 06 07 08 05 
S 


6 差 流 波形 和 判 据 HS 值 计算 结果 〈 算 例 4) 
Fig.6 Differential current and value of HS for the proposed 


criterion(Scenario 4) 


图 7 差 流 二 次 谐 波 占 基 波 百分比 〈 算 例 4) 


Fig.7 The ratio of second harmonic to fundamental of 


differential current (Scenario 4) 


根据 图 6， 带 高 阻 故障 空 载 合 闸 情 况 下 ， 在 判 
据 启 动 后 ， 涌 流 特征 明显 的 正 半 周 ，Zs 值 低 于 门 
槛 值 0.77， 保 护 被 短暂 闭锁 ;但 到 涌流 值 较 小 的 
负 半 周 ，HS 值 迅速 增 大 并 超过 0.77 的 动作 门槛 值 ; 
可 以 看 到 在 人 0.4105s， 即 变压器 带 高 阻 故障 合 闸 
后 约 半 个 周波 后 ， 判 据 即 解除 闭锁 ， 开 放 保 护 ， 
正确 动作 。 对 于 此 类 差 流 为 故障 电流 县 加 涌流 的 情 
况 ， 差 流 二 次 谐 波 百分比 可 能 会 高 于 15% 的 制 动 
门槛 (如 图 77， 二 次 谐 波 制 动 判 据 将 误 闭锁 保护 
《或 造成 保护 延 时 动作 ) ; 但 采 
则 能 在 很 短 时 间 内 开放 保护 ， 使 保护 正确 动作 。 

限于 篇 幅 ， 以 下 算 例 5-7 只 给 出 差 流 波形 和 
Hausdorff 距离 判 据 HS 计算 值 序列 。 

算 例 5: 该 算 例 为 空 载 合 闻 励 磁 涌 流 后 又 发 生 
变压器 内 部 故障 的 场景 ，{=0.4s 变压器 空 载 合 闸 ， 
{=0.605s 发 生 Y/Y 变 一 次 侧 出 口 三 相 接 地 故障 。 


j 本 文 所 提 判 据 ， 
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8 差 流 波形 和 判 据 ZS 值 计 算 结果 〈 算 例 5) 
Fig.8 Differential current and value of HS for the proposed 


criterion(Scenario $) 


如 图 8 可见 ， 在 空 载 合 闻 典 型 励磁 涌流 阶段 ， 


HS 计算 值 一 直 稳 定 地 低 于 0.77 的 动作 门槛 ， 保 护 
被 可 靠 闭锁 ; 当 内 部 故障 发 生 后 ， 在 f=0.61s 时 刻 


(内 部 故障 发 生 约 /4 个 周波 ) HS 计算 值 超过 
0.77 的 动作 门槛 ， 解 除 闭 锁 ， 保 护 立刻 正确 动作 。 
算 例 6: 该 算 例 为 变压器 经 历 区 外 故障 转换 为 
区 内 故障 的 场景 ， 伍 0.305s 发 生 Y/Y 变 外 部 三 相 短 
路 故障 ， 人 0.505s 故障 转 为 Y/Y 变 一 次 侧 出 口 三 相 
接地 故障 。 
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9 差 流 波形 和 判 据 25 值 计算 结果 〈 算 例 6) 
Fig.9 Differential current and value of HS for the proposed 


criterion(Scenario 6) 


区 外 故障 阶段 ， 差 流 幅 值 较 小 ， 判 据 不 予 启 
动 ; 当 区 外 故障 转 为 变压器 区 内 故障 后 ， 判 据 立 
刻 启 动 ， 在 {=0.51s， 即 故障 转 为 区 内 故障 后 1/4 个 
周波 ，HS 计算 值 超过 动作 门槛 ， 保 护 迅 速 正确 动 
作 。 

从 上 述 仿真 算 例 的 结果 可 以 看 到 ， 本 文 所 提 
基于 Hausdorff 距离 算法 的 差 动 保护 新 判 据 ， 在 应 
对 变压器 经 历 各 种 内 部 故障 、 励 磁 涌流 和 发 展 性 故 
障 扰动 时 ， 均 能 够 做 出 正确 判断 ， 并 在 速 动 性 方 
面 优势 明显 。 

算 例 7: 抗 干扰 性 分 析 


根据 前 文 分 析 ，Hausdor 任 距离 算法 主要 基于 于 Hausdorff 距 离 算 法 新 判 据 在 抗 干 扰 性 方面 有 较 
序列 的 整体 特征 进行 判断 ， 因 此 ， 受 差 流 中 非 周 。” 强 优 势 。 
期 分 量 和 高 次 谐 波 的 影响 很 小 ， 并 且 在 整体 特征 。 7 结论 
相同 的 情况 下 部 分 点 缺失 对 Hausdorff 距离 算法 的 
影响 较 小 ， 因 此 ， 该 算法 具备 抗 数据 丢失 的 性 能 ， 本 文 基于 涌流 和 故障 电流 波形 整体 形态 的 特 
可 利用 该 性 能 ， 无 需 增 加 任何 滤波 环节 即 可 使 判 。” 征 差 异 ， 结 合 Hausdorff 距 离 算 法 在 波形 相似 性 判 
据 具备 非常 好 的 抗 干扰 能 力 。 具 体 实现 方法 是 : 在 。” 列 中 的 优势 ， 提 出 了 基于 Hausdorff 距 离 算法 的 变 
进行 Hausdorff 距离 计算 前 ， 根 据 实际 情况 设计 在 。 压 器 关 动 保 扩 新 判 据 ， 并 进行 了 仿真 试验 验证 。 结 
每 个 时 间 窗 内 舍弃 若干 极 值 点 一 一 对 于 没有 噪声 。” 末 表 明 ， 本 文 所 提 新 判 据 能 够 快速 正确 地 区 分 变 
干扰 的 差 流 序列 ， 使 弃 部 分 点 不 影响 整体 波形 特 ”上 压 器 励磁 涌流 和 故障 差 流 ， 对 于 二 次 谐 波 和 波形 
征 ， 因 此 对 Hausdo 人 ff 距离 计算 值 几乎 没有 影响 ; 对 称 判 据 容易 误 判 的 对 称 性 涌流 ， 也 有 较 强 的 识 
对 于 有 噪声 干扰 的 差 流 序列 ， 若 每 个 时 间 窗 内 含 ”， 别 能 力 ， 可 防止 差 动 保护 误 动 ， 并 能 有 效 识别 故 
弃 1 个 极 值 点 ， 则 判 据 在 采用 1/4 周波 时 间 窗 时 ， 障 电 流 谷 加 涌流 的 情况 。 该 算法 的 优势 还 主要 体现 
一 个 周期 差 流 序列 即 可 识别 并 舍弃 至 少 4 个 噪声 干 。 “在 两 个 方面 ， 其 一 是 计算 量 小 ， 判 别 迅速 ， 不 会 
扰 点 ， 足 够 满足 工程 实际 中 一 个 周波 内 最 多 考虑 。 ”对 保护 启动 造成 其 他 延 时 ， 最 快 能 达到 跟 保 护 启 
一 个 噪声 干扰 点 的 要 求 。 动 同 时 判别 ， 一 般 故 障 情 况 5ms 便 能 做 出 判断 ; 
二 二 后 丽 区 放流 其 二 因 Hausdorff 距 离 算法 抗 数据 丢失 能 力 强 ， 因 


i 此 对 于 差 流 波形 中 有 个 别 噪声 干扰 的 情况 ， 在 执 
(ol 行 算法 具体 计算 前 ， 设 置 舍弃 部 分 极 值 点 ， 不 影 
| 响 判别 结果 ， 体 现 了 该 算法 优越 的 抗 干扰 性 能 。 


T 


由 于 励磁 涌流 波形 特征 受 剩 磁 、 绕 组 连接 方式 
人 合 阅 相 角 等 因素 的 影响 ， 以 波形 特征 为 依据 构成 


的 判 据 及 其 门槛 值 的 选取 都 无 法 保证 保护 在 所 有 
情况 下 都 不 误 动 ， 就 如 二 次 谐 波 制 动 判 据 一 样 ， 


0 0.35 0.4 045 05 0.55 0.6 随 着 技术 进步 和 制造 工艺 的 提升 ， 即 使 将 二 次 谐 
EE 波 百分比 门槛 设置 在 一 个 很 低 的 值 ， 也 无 法 确保 
和 | 所 有 涌流 扰动 时 差 动 保护 都 能 可 靠 制 动 。 本 文 所 提 
: 判 据 ， 是 对 基于 电流 波形 特征 应 用 的 新 的 尝试 ， 
| 为 变压器 差 动 保护 算法 提供 了 一 种 新 的 实现 方式 ， 
办 后 续 将 对 算法 有 效 性 进行 动 模 试验 和 现场 数据 的 
图 10 差 流 波形 和 判 据 HS 值 计算 结果 〈 算 例 7) 验证 ， 并 借助 大 量 数据 选择 合理 的 整定 值 ， 不 断 
Fig.10 Differential current and value of HS for the proposed 对 算法 进行 优化 和 完善 。 
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